Two-chamber lattice model for thermodiffusion in polymer solutions 
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When a temperature gradient is applied to a polymer solution, the polymer typically migrates 
to the colder regions of the fluid as a result of thermal diffusion (Soret effect). However, in recent 
thermodiffusion experiments on poly(ethylene-oxide) (PEO) in a mixed ethanol/water solvent it is 
observed that for some solvent compositions the polymer migrates to the cold side, while for other 
compositions it migrates to the warm side. In order to understand this behavior, we have developed 
a two-chamber lattice model approach to investigate thermodiffusion in dilute polymer solutions. 
For a short polymer chain in an incompressible, one-component solvent we obtain exact results for 
the partitioning of the polymer between a warm and a cold chamber. In order to describe mixtures 
of PEO, ethanol, and water, we have extended this simple model to account for compressibility 
and hydrogen bonding between PEO and water molecules. For this complex system, we obtain 
approximate results for the composition in the warmer and cooler chambers that allow us to calculate 
Soret coefficients for given temperature, pressure, and solvent composition. The sign of the Soret 
coefficient is found to change from negative (polymer enriched in warmer region) to positive (polymer 
enriched in cooler region) as the water content of the solution is increased, in agreement with 
experimental data. We also investigate the temperature dependence of the Soret effect and find 
that a change in temperature can induce a change in the sign of the Soret coefficient. We note a 
close relationship between the solvent quality and the partitioning of the polymer between the two 
chambers, which may explain why negative Soret coefficients for polymers are so rarely observed. 



I. INTRODUCTION 

A temperature gradient applied to a fluid mixture gen- 
erally induces a net mass flow, which results in the for- 
mation of a concentration gradient. This effect is known 
as thermodiffusion or the Ludwig-Soret effect. 1 ! 2 ! 3 ! 4 The 
Soret coefficient St relates the steady state concentration 
gradient to the imposed temperature gradient. By con- 
vention, the Soret coefficient of component i is positive 
if component i is enriched in the cooler region. 5 Thermal 
diffusion has long been used as an effective tool for sep- 
arating mixtures of isotopes^ More recently, the effect 
has been used to char acter ize mixtures of complex fluids 
(see for example Refs. 



Thermodiffusion in a binary fluid mixture is described 
by the flux of one of the components in response to a tem- 
perature and concentration gradient i The flux is given 
by 



Ji = -pDVcx -pci(l-ci)X>'Vr, 



(1) 
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In liquid mixtures whose components differ widely in 
molecular mass, such as polymer solutions^ and col- 
loidal suspensions^, it is typically the heavier component 
that migrates to the cold region. There are, however, 
exceptions. In 1977, Giglio and Vendramini found a neg- 
ative Soret coefficient for poly(vinyl alcohol) in water^ 
Very recently, de Gans et ali 10 ^ 11 reported results of ther- 
mal diffusion forced Raleigh scattering (TDFRS) mea- 
surements on solutions of poly(ethylcne oxide) (PEO) 
in mixtures of ethanol and water. In pure water, PEO 
shows the expected migration to the cold region of the 
fluid (St > 0). However, in solutions with low water 
content, PEO is found to migrate to the warmer region 
of the fluid (St < 0). Although changes in sign of the 
Soret coefficient have been reported for a number of liq- 
uid mixtures of small- molecule fluids, including alcohol 
solutions jl^iiiiiii 5 . the PEO/ethanol/ water system ap- 
pears to be the first polymer solution for which such a 
sign change has been observed. 



where D is the mutual diffusion coefficient, D' the ther- 
mal diffusion coefficient of component 1, p the total mass 
density, c\ the mass fraction of component 1, and T is 
the temperature. Here the pressure is assumed to be con- 
stant throughout the mixture and the flux J± describes 
the flow of component 1 with respect to the center of 
mass of the system. 1 Eventually, the system reaches a 
stationary state in which the flux J\ vanishes. Inserting 
Ji = into Eq. Q yields 



1 Vci D' 



ci(l-ci)VT D 



(2) 



The Soret coefficient of component 1 is the ratio of ther- 
mal and mutual diffusion coefficients 



s - D ' 



(3) 



More generally, we define the Soret coefficient of compo- 
nent i of a mixture as 



1 



Ci(l - Ci) dT 



(4) 



In order to describe thermodiffusion in ternary mix- 
tures, such as PEO in a mixed solvent, concentra- 
tion gradients and fluxes of two of the components are 
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considered 

Ji = -p£) u Vci- J o£> 12 Vc2-pci(l-ci) J D / 1 VT, (5) 
J% = -pD 2 iVa- pD 22 Vc 2 - pc 2 (l-c 2 )D' 2 VT, (6) 

where J\ and J2 are the flux of component 1 and 2, re- 
spectively, and the flux of the third component is given 
by J3 = — Ji — J 2 . Due to Onsager's relations, only three 
of the four isothermal diffusion coefficients Dij are in- 
dependent. Their relationship has a complicated form, 
however, since the fluxes are defined with respect to the 
barycentric reference framei. The prefactors of the gen- 
eralized thermal diffusion coefficients D\ have been cho- 
sen in such a way that Eq. (JJJ for the binary mixture is 
recovered for c 2 — or ci = 0. In the steady state, each 
of the fluxes vanishes and the relationship @ between 
composition and temperature gradient generalizes to 
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(7) 



and similarly for component 2. For ternary mixtures, we 
continue to define the Soret coefficient of component i 
through Eq. J3J . 

Thermal diffusion in liquid mixtures is not well under- 
stood and even the sign of the Soret coefficient cannot 
generally be predicted (see e.g. Refs. 000)- Due to the 
complexity of the task, attempts to extend the kinetic 
gas theoryii of thermodiffusion to the liquid state have so 
far been unsuccessful 4 *^. Molecular dynamics simulations 
(for a review see Ref. ITsf) have become an important tool 
in the investigation of thermodiffusion in small-molecule 
liquids. Long computation times make it difficult, how- 
ever, to address thermodiffusion in polymeric systems. 

Earlier theoretical work on thermodiffusion, see e.g. 
Refs. llll2TlllTMTl2^ll2lll2ll2^. made use of the "heat of 
transport" concept. For dilute binary solutions, the Soret 
coefficient of the solute is related to the heats of transport 
of the two species as follows^ 3 - 



St — 



Ql Q 2 
fc R T 2 : 



(8) 



where is Boltzmann's constant. The numerator rep- 
resents the net energy flowing through the boundary of 
a small volume of the mixture when the composition 
inside the region changes while temperature and pres- 
sure remain constant. Various microscopic expressions 
for the heat of transport have been obtained with the 
aid of lattice model calculations. Wirtz and HibyiliSfliSi 
considered two contributions to the energy required for 
a particle to move from one location to another, the 
energy to detach the particle from its neighbors and 
the energy to create a void for the particle to move 
into. Denbigh 23 considered the energy associated with 
particle diffusion across a boundary to be a "certain 
fraction" of the change in potential energy, which is 



represented by nearest-neighbor interactions in random- 
mixing approximation^ 3 - This type of approach is use- 
ful in describing thermal diffusion in regular solutions 4 
but fails for more complex systems. Prigogine et alM- 
found that a model based on energies of nearest-neighbor 
interactions alone^ led to contradictions when applied 
to alcohol mixtures, where the Soret coefficient changes 
signiSii 3 .. They argued that a "free energy" for detach- 
ing a molecule needs to be considered in systems with 
associating solvents. By assuming that alcohol molecules 
form complexes whose disruption results in a loss of local 
energy and a gain in local entropy, they were able to ex- 
plain qualitatively the sign change of the Soret coefficient 
for these systems^ 

In this work, we investigate the Soret effect in di- 
lute polymer solutions with the aid of two-chamber 
lattice models. Following traditional experimental 
methods^Si&l we consider a system divided into two 
chambers of equal size that are maintained at slightly 
different temperatures. Particles arc free to move be- 
tween the chambers, which do not otherwise interact. If 
the pressure differences between the chambers are small 
enough to be neglected, the Soret coefficient can be deter- 
mined from the difference in composition of the solutions 
in the two chambers ^iVi 4 We start by investigating a 
polymer chain in an incompressible, one-component sol- 
vent in Section [n] This allows us to introduce the gen- 
eral strategy and obtain exact results for a simple model 
system. Section lnTl focuses on PEO/ethanol/ water mix- 
tures. In Section IlII Al we present a lattice model for the 
static properties that includes the effects of compress- 
ibility and specific interactions between PEO and water. 
The corresponding two-chamber system is introduced in 
Section UTTBI Results of our calculations are presented in 
Section llVl and compared with experimental data, where 
available. In Section [V] we discuss the work presented 
here and provide an outlook to future work. The ap- 
pendices describe the determination of system-dependent 
parameters for PEO/ethanol/ water mixtures and the re- 
lation between physical properties and model variables. 



II. LATTICE MODEL FOR INCOMPRESSIBLE, 
DILUTE POLYMER SOLUTIONS 



A. Single chamber 

We consider a simple cubic lattice of N sites. Since we 
are interested in dilute solutions, we assume that only 
one polymer chain is present. A chain of N c beads has 
N c — 1 bonds and occupies without overlap N c contigu- 
ous lattice sites. For an incompressible system, all of 
the remaining lattice sites are occupied individually and 
represent solvent molecules. The filling fraction of the 
polymer and the solvent are defined as 4> p — N c /N and 
4> 8 = N s /N, respectively, with <p s = 1 — <j) p . Interactions 
between non-bonded chain segments on nearest-neighbor 
lattice sites are described by an interaction energy e pp . 
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net repulsive polymer-polymer interactions 

Ae = +1 




net attractive polymer-polymer interactions 

Ae = -1 



10 



FIG. 1: Radius of gyration squared of chains of length 
iV c = 17 as a function of reduced temperature for net attrac- 
tive (Ae < 0) and net repulsive (Ae > 0) polymer-polymer 
interactions. 



r 


1 { 


b 






\ i 


) 




r 










I i 


s 




J \ 


r 






J ^ 






J % 


i ' 




f — % 


t f 
■> r 


J 


( 


) — c 


) — < 
















£ pp' 






J L 


) 












) — c 


) 

i 


c 


) ( 

In, 




►— i 


> — < 


HH 


► 


D 


► — < 




^ r 






i — < 


> — ( 


► — < 


1 — < 


) 

3 










\j \-j » 


►— < 


>— < 


> — C 



FIG. 2: Two-dimensional illustration of the two-chamber 
lattice model for the incompressible case. The right chamber 
contains the polymer chain, where the black circles and heavy 
lines indicate the beads and bonds of the polymer chain while 
the grey circles represent solvent sites. 



each conformation as well as the radius of gyration. The 
radius of gyration for a given temperature is calculated 
from 24-25 



Polymer-solvent and solvent-solvent interactions are also 
restricted to nearest-neighbor lattice sites and described 
by interaction energies e ps and e ss , respectively. 

The number of chain conformations is given by the 
number of self-avoiding random walks of N c — 1 steps 
and can be enumerated exactly for short chains (see e.g. 
Refs. I24l2qh The probability for a given chain confor- 
mation to be realized depends on the interactions in the 
system. For a chain conformation with m pair contacts 
between polymer segments, the internal energy of the 
polymer-solvent system is given by 
' z 
v2 



E(m) 



I) Nr. 



lj Ae 

((f-DiV. + l.e 



pp 



(9) 

where z = 6 is the coordination number of the simple 
cubic lattice and Ae = e pp + e ss — 2e ps is the net inter- 
action energy between polymer segments. The canonical 
partition function of the system is given by 



Zp i(T) =NJ2 c(m) exp(-/3£(m)), 



(10) 



where c(m) is the number of chain conformations with 
m polymer-polymer contacts and (3 = 1/ksT, where T 
is the temperature and fee is Boltzmann's constant. The 
internal energy for a given temperature is calculated in 
the usual way 



C^ P ol 



iV^c(TO)£ , (m)exp(-/3£;(m))/Z po i (11 

rn 

.2<9m(Z P oi) 



dT 



(12) 



In a preliminary step, we generated all single chain 
conformations of chains up to length N c — 17 and de- 
termined the number m of polymer-polymer contacts for 



N c(m)R2(m) exp( 



-f3E(m))/Z poh (13) 



where i? 2 (m) is the average radius of gyration squared 
for conformations with m contacts. For the incompress- 
ible system, we express all interaction energies in terms 
of a positive energy unit e, which also sets the tem- 
perature scale; the reduced temperature T* is given by 
T* = k^T/e. In Fig. ^ we present chain dimensions for 
chains of length iV c = 17 as a function of reduced tem- 
perature for the case of net attraction (Ae < 0) and net 
repulsion (Ae > 0) between the polymer segments. As 
expected, the chain expands with increasing temperature 
for the case Ae < 0, while the chain decreases in size with 
increasing temperature for the case Ae > 0. 

For polymers in dilute solutions, the chain dimen- 
sions are an indicator for solvent quality^; the better 
the solvent the larger the chain dimensions. For very 
long chains with net attractive interactions between the 
polymer segments, a collapse transition is observed at 
the 8 temperature of the polymer solution^. Monte 
Carlo simulations 2 ^ yield an estimate of 9* = 3.60 ± 0.05 
for the reduced 9 temperature of an infinitely long iso- 
lated chain. From our exact enumeration results we de- 
termined the dimensions of an isolated chain of length 
N r . = 17 at the 9 temperature of the infinite chain and 
) =4.13 in units of lattice constant squared. 



found Rg 



B. Two-chamber system 

In order to investigate thcrmodiffusion, consider a 
single-chain system that is divided into two chambers, 
A and B, with slightly different temperatures, T A > T B , 
see Fig. [21 Under the assumption that the chambers are 
non-interacting, the partition function of the whole sys- 
tem is the product of the partition functions of the indi- 
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vidual chambers, Z A Z B . The chambers are represented 
by lattices with N a , a <E {A, B} sites. The canonical par- 
tition function of the chamber that contains the polymer 
is given by Eq. (|1(J|I . the canonical partition function for 
the chamber without polymer consists of a single term 



J nop 



(T ) = exp(-/3 a £ , nop ) 



(14) 

The internal 



with a e {A,B} and E nop = §A a e SE 
energy in this case is independent of temperature, U nop 



nop ■ 



If the polymer chain is allowed to be in either 



chamber, the sum of states for the system is given by 



Q = Zp \(TA)Z nop (TB) + Zp i(TB)Z nop (TA 



(15) 



Accordingly, the probability q A to find the polymer in 
the (warmer) chamber A can be calculated from 



<?a = Zp \(TA)Z n0 p(TB,)/Q. 



(16) 



qB = 1 — <za is the probability to find the chain in the 
cooler chamber B. 

For small temperature differences ST = Ta — Tb we set 
T B = T,T A = T + 5T and rewrite Eq. HH as 



1 



= 1 



Z pol (T)Z nop (T + ST) 



<7A Zp oi (T + ST)Z n0 p(T) ' 

Expanding the partition functions according to 



(17) 



ln(Z(T + ST)) = ln(Z(T)) + f3U(T)ST/T, (18) 



where (3 — 1/fceT, we find for the probability 



1 / , x ST 

= 1 + CXp j3 (U n0 p ~ Upol) — 

QA V 1 



(19) 



We employ temperature differences ST JT = 10~ 4 in our 
calculations and find the results from Eqs. lfTfi|l and l(T§|l 
to be indistinguishable. In order to compare with the the- 
ories based on the heat of transport) 1 ! 2 ! 3 ! 4 ! 13 ! 19 ' 20 ! 21 ! 22 ! 23 
we expand the expression for the probability qA — 1/2 
from Eq. I|19|) in powers of ST /T and obtain to first or- 
der 

gA-~^-~/^nop-£W^. (20) 

For this simple model, the probability to find the polymer 
in the warmer chamber A is larger than 1/2 when the in- 
ternal energy of the chamber without polymer is smaller 
than the internal energy of the chamber with polymer at 
temperature T . Hence, the difference in internal energy 
between two chambers, one with and one without poly- 
mer, at the same temperature determines the probability 
to find the polymer in the warmer of two chambers. 

The Soret coefficient ( Eq. (0J ) can be calculated from 
the filling fraction as well as from the mass fraction, since 
dci/dipi = c\Cij '4>\4>2- Setting 4>i = p , we obtain in the 
dilute limit 
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OpA 



9 pBj 



(T A - T B ) 



(21) 
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FIG. 3: Soret coefficients of an incompressible dilute poly- 
mer solution as a function of temperature. The curves repre- 
sent exact enumeration results for chains of length N c = 17, 
where the Soret coefficient has been calculated from Eq. 1231 . 
The interaction parameters are as indicated in the figure and 
correspond to net repulsive interactions (Ae = +1) for the 
dashed curve and net attractive interactions (Ae = — 1) for 
the dotted and solid curves, respectively. 



where the filling fractions in the chambers are calculated 
from the probability <7a 

P A = qAN c /N A , pB = (1 - q A )N c /N B . (22) 

Typically we assume that the chambers are equal in size 
and set Aa = Ab = N/2, which yields for the Soret 
coefficient 
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,«A 



ST 



(23) 



With the aid of Eq. (|2U|) for the probability, the Soret 
coefficient becomes 



St 



U n 



Upol 



k B T 2 



(24) 



In agreement with the definition of the Soret coefficient, 
this expression for St is independent of the tempera- 
ture difference ST. It is also independent of the poly- 
mer concentration, in agreement with ex per imental re- 
sults for dilute polymer solutions (cf. Ref . l28i) . Compar- 
ing Eq. H24fl with the heat of transport expression for 
dilute solutions in Eq. ©, we see that the internal en- 
ergy difference U nop — U po \ takes the place of the heat of 
transport difference Q\ — Q* 2 . 

Results for the Soret coefficient calculated from 
Eqs. {THJl and are presented in Fig. We note 
first that both positive and negative Soret coefficients 
are realized and that St may change sign as a func- 
tion of temperature. For the simple model discussed 
in this section, the internal energy £/ nop of a solvent- 
filled lattice is independent of temperature. If U po i was 
also independent of temperature, then St would vary as 
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l/T 2 , according to Eq. l|2"Ifl. The deviation of the ac- 
tual calculated Soret coefficients from this l/T 2 depen- 
dence, most notably the sign change of St, is due to the 
temperature variation of the internal energy U po i of the 
polymer system. A random mixing approximation for 
polymer-polymer and polymer-solvent contacts yields a 
temperature-independent internal energy U po \. Hence, in 
order to calculate Soret coefficients of polymer solutions, 
one has to go beyond the random mixing approximation, 
as we do in this work. 

A comparison of the results for the chain dimensions 
in Fig. H with the results for the Soret coefficients rep- 
resented by the dashed and solid curves in Fig. shows 
that the Soret coefficient decreases as the solvent quality 
increases. We find this trend confirmed when we consider 
the more complex system below. However, the Soret co- 
efficients represented by the dotted curve show that this 
need not be the case. The chain dimensions for this case 
are the same as for the case represented by the solid line, 
since they correspond to the same net interaction energy 
Ae = — 1. It was recognized early 22,23 that, in contrast to 
thermodynamic and configurational properties, thermod- 
iffusion is not determined by the net interaction energy 
Ae but by all three interaction parameters. In fact, the 
simple heat of transfer model a 21 i 22 i 23 are most successful 
for regular solutions, where Ae = 0. 



III. LATTICE MODEL FOR PEO IN 
ETHANOL/ WATER MIXTURES 

Solutions of high molecular weight poly (ethylene ox- 
ide) (PEO) in ethanol and water have interesting prop- 
erties. Hydrogen bonding between PEO and water 
molecules plays an important role in aqueous solutions 
of PEO (see e.g. Refs. l29ll3Cf) . Water is a good solvent 
for PEO at standard temperature and pressure. How- 
ever, the solvent quality decreases with temperature and 
a miscibility gap opens above a lower critical solution 
temperature^ Ethanol, on the other hand, is a poor 
solvent for PEO at room temperature but the solubility 
increases with temperature^. In mixtures of ethanol and 
water at standard temperature and pressure, the water 
content determines the solubility of PEO. For the molecu- 
lar weight considered in this work, the transition between 
poor and good solvent condition appears between a water 
content of 5% and 10% by weighti2iii. Light scattering 
experiments- show that the PEO chains expand with 
increasing water content, indicating that the addition of 
water improves the solvent quality. 



A. Static properties 

In order to describe dilute solutions of poly(ethylene 
oxide) (PEO) in mixtures of ethanol and water, we ex- 
tend the lattice model introduced above. At a given tem- 
perature, pressure, and composition, the single chamber 



system is represented by a simple cubic lattice with N 
sites, of which 7V C , TVs, and _/V w are occupied by the poly- 
mer (PEO), the first solvent (ethanol), and the second 
solvent (water), respectively. In order to account for 
compressibility, we allow sites to be unoccupied so that 
N = N c + N s + N w + N v , where 7V V is the number of 
voids. The total volume of the lattice is V — vqN, where 
vq is the volume of one elementary cube. 

Interactions between occupied nearest neighbor sites 
are described by interaction energies e^ , where the sub- 
scripts indicate the occupants of the sites (p for poly- 
mer, s and w for the solvents; voids are assumed to have 
zero interaction energies). In aqueous solutions, hydro- 
gen bonding between PEO and water plays an important 
role (cf . Ref . 12^) . In order to account for these specific 
interactions, we introduce an orientational degree of free- 
dom in the description of water. Each elementary cube 
occupied by water is assumed to have one special face. 
If this face is exposed to a polymer segment, the interac- 
tion energy is e pw;s (strongly attractive) otherwise e pw;n 
(non-specific). 

Because of the complexity of the system, we can no 
longer evaluate the canonical partition function exactly. 
In the following, we employ the exact enumeration re- 
sults for the chains of length iV c = 17 and consider ex- 
plicitly contacts that involve the polymer but we em- 
ploy a random mixing approximation for all other con- 
tacts. A chain conformation with m pair contacts has 
n n = 4iV c + 2 — 2m nearest neighbor (nn) sites, which are 
occupied by n^, i € {s,w, v} solvent particles and voids. 
With the aid of the random mixing approximation for all 
but the polymer contacts, the canonical partition func- 
tion of the system can be written as 



Z pol (N,T,N Wl N s ) 



(25) 



JVXc(m)£6^-"* 



N — n n — N c 
N w - n w 



[n 8 ] 



N — n n — N c — (JV W - n w ) 
N s - n R 



, e -/3(me pp +n s e ps ) ^g e -/3e pWi n _|_ e -/3e pw;s y 



-PEt 



where, as before, c(m) is the number of chain confor- 
mations with m polymer-polymer contacts. The square 
brackets around the summation indices indicate that the 
summation is performed consistent with the available 
nearest neighbor sites and the total filling of the lattice. 
The energy E x denotes the contribution to the total en- 
ergy due to solvent-solvent interactions evaluated in ran- 
dom mixing approximation^, cf. Eq. (|29|l below. 
The pressure of the system is calculated from 



P = 



k B T f d In Zppi 
dN 



vo 



(26) 



N B ,N„,N C 



By performing partial summations over the terms in 
Eq. (|25|l the probabilities for specific sets of states can 
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FIG. 4: Radius of gyration squared, Rg, as calculated from 
Eq. I28H . Panel (a) shows chain dimensions as a function 
of solvent composition at temperature T — 293 K, pressure 
P « 0.1 MPa, and a PEO concentration of 5 g/L. Panel (b) 
shows the temperature variation of the chain dimensions when 
the solutions are heated at constant pressure. The dashed line 
indicates the chain dimensions, R g (8*), of the isolated chain 
at the 9 temperature of the infinite chain, as discussed in 
Section lira 



be determined. If we write the partition function as 

Z p ol — ^ ^m,n w ,n s j (^^) 

m,[n w ] , [n s ] 

the average radius of gyration (Rg) is given by 

( R D = Z poI E Rl(m)Z m , nw<ne . (28) 



In Appendix we discuss how the system-dependent 
parameters of the lattice model were determined. In Ap- 
pendix [B] we describe how lattice size and occupation 
numbers for given temperature, pressure and composi- 



tion of the mixture are calculated. Since chain dimen- 
sions are an indicator for solvent quality, we employed 
calculated Rg values to obtain estimates for the mixed in- 
teraction parameters. In Fig. ^ we present graphs for the 
chain dimensions calculated with the aid of the system- 
dependent parameters presented in Table |U 



B. Two-chamber system 



As in section III Bl we consider a system of two cham- 
bers with slightly different temperatures and determine 
the probability to find the polymer in the warmer of the 
two chambers. Assuming again that the chambers are 
non-interacting, the canonical partition function of the 
system for a given occupation of the chambers is the 
product of the individual partition functions. 

To find the required canonical partition function of a 
chamber without polymer, consider a lattice of N sites 
occupied by the two types of solvent and voids, N = 
N s + N w + N v . Denoting the filling fractions of ethanol 
and water by ^> s = N s /N and ^> w = A w /A and assuming 
random mixing, the internal energy is given by2i 



E, 



"•>p = ^ N ( £ss< ^ + £v 



+ 2e 



wsrsrw I 



(29) 



Accordingly, the canonical partition function of the lat- 
tice without polymer takes the form 



Z nop (N,T,N s ,N w ) 



G 



N- A w 



(30) 



To ease notation for the two-chamber sum of states, 
we define the canonical partition for a single chamber as 



Z(N,T,N w ,N s ,N p ) 

' Z pol (N,T : N s ,N w ) for N p = 1 
Z nop {N 7 T,N s ,N w ) for N p = 



(31) 



where N p € {0, 1} is the number of polymer chains in the 
chamber. 

For the mixed solvent system considered here, the con- 
centration of the two solvents will generally be different 
in the warm and cold regions of a fluid. Since we have no 
a priori information about concentrations (or chemical 
potentials) of the solvents, we consider all distributions 
of particles consistent with fixed total particle numbers. 
The sum of states is then given by 



Q = E E E 

JV p , A =0 [AT W _ A ] [Af s-A ] 

x Z(N-N A ,T B ,N, 



Z(N A ,T A ,N WtA ,N B , A ,N p , A ) 
, -N W . A ,N S - A S , A ,1 



(32) 



where, as before, square brackets indicate summations 
consistent with the total numbers of particles and lat- 
tice sites. Also as before, chamber A is considered the 
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warmer chamber so that ST = T A — Tb > 0, and equal- 
sized chambers are used, N A — N% = N/2. As we are 
performing the calculation of the terms in the sum of 
states, we monitor for each chamber the composition and 
the pressure of the mixtures. This allows us to calculate 
the average quantities for each chamber by performing 
weighted sums. For example, the average mass fraction 
of component i in chamber A, is calculated from 



l 1 

C *. A = n E E E 

^ AT PiA =o [jv w ,a] [JV s ,a" 



Ci(N w , A ,N s ,A,N p , A ) (33) 



X Z(Ar A ,T A ,iV w ,A,^ s> A,iV P; A) 
x Z{N - N A , T B ,N W - JV Wj a, N s 



A S , A , 1 - N pA ), 



where Cj is the mass fraction of component i calculated 
according to Eq. (|B4(1 . 

The probability to find the polymer in chamber A in 
this model is given by 



<1A 



\ E E ^poi(A a ,T a ,A w ,a,A S: 



Q 



(34) 



[JVw.a] [N.,a] 



x Z nop (N - N A , T B ,N W - A WiA , N s - N S<A ). 

As in the case of the incompressible system, this prob- 
ability is related to an internal energy difference of two 
chambers at the same temperature: 



<7A 



1 (t/nop) - (Upol) ST 

'4 



k B T 



T 



(35) 



The angular brackets indicate an average over all configu- 
rations of particles in two chambers at the same tempera- 
ture T, where the polymer is confined to one of the cham- 
bers. C/p i and £/nop are the average internal energies of 
the chamber with and without polymer at fixed composi- 
tion. For the small temperature differences ST = 1CP 4 K, 
corresponding to ST /T ~ 3 x 10 -7 , employed in our cal- 
culations, values for the probability q A calculated from 
Eqs and (|33|) agree to more than five digits. As in 
the incompressible case, the excess probability q A — 1/2 is 
proportional to the temperature difference and indepen- 
dent of the polymer concentration for dilute solutions. 
Relations similar to Eq. I1H5II can be derived for the prob- 
ability of finding a given number of solvent particles in 
the warmer chamber A. They show that the difference 
between the average numbers of particles in chambers A 
and B is proportional to ST. 



IV. RESULTS 

We have applied our lattice model to solutions of PEO 
in ethanol/ water mixtures under a variety of conditions. 
As noted above, calculated values for composition dif- 
ferences in the two chambers are expected to be propor- 
tional to the temperature difference ST = T A — T B for suf- 
ficiently small values of ST. We find this to be the case 
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FIG. 5: Excess probability to find the PEO chain in the 
warmer chamber A as a function of water content of the so- 
lution. Calculations were performed for the compositions in- 
dicated by the symbols; the symbol size is larger than the 
uncertainty due to discrete occupation number discussed in 
Appendix IbI 
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FIG. 6: Soret coefficient of PEO in mixtures of ethanol and 
water at room temperature and atmospheric pressure. The 
filled symbols represent experimental data by de Gans et alw—, 
the open symbols represent preliminary experimental results 
by Kita et aim— , the line is the result of our lattice model 
calculations. 



for a large range of ST values and used ST = 10~ 4 K 
throughout the calculations. Furthermore, since the 
probability to find the polymer in the warmer chamber 
is independent of polymer concentration, we may choose 
the PEO content of the solution. As discussed in Ap- 
pendix [51 a low polymer concentration corresponds to 
a large lattice size, which reduces the uncertainty due 
to integer occupation numbers but increases the compu- 
tation time. We found that the PEO concentration of 
5g/L used in the experiments represents a satisfactory 
compromise. The results presented in Figs. pertain 
to room temperature (293 K), atmospheric pressure (0.1 
MPa), and a PEO concentration of 5 g/L. 

In Fig. [S] we present values for the excess probabil- 
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FIG. 7: Composition differences Set = Cj,A — c^.b between 
chambers A and B for the three components (PEO, ethanol, 
water) as a function of water content of the solution. The 
compositions of the chambers were calculated from Eq. I.SMI 
and divided by the temperature difference ST = Ta — Tb = 
1(T 4 K. 



ity <7a — 1/2 as calculated from Eq. (|34|l . In Fig. we 
present the corresponding values for the Soret coefficient 
of PEO calculated according to Eqs. ffify and Q. For 
comparison, we include experimental data by Wiegand 
and coworkersifl^. Both experiment and theory show a 
change in sign of the Soret coefficient as the water con- 
tent of the solution is increased. For low water concen- 
trations of the solution, the polymer is more likely to be 
found in the higher temperature chamber; for high water 
concentrations the opposite is true. Differences between 
theory and experiments are most pronounced at low wa- 
ter concentrations, where our calculations overestimate 
the Soret effect. This is a consequence of our choosing a 
mixed interaction parameter e ps that emulates for short 
chains the poor solvent conditions that long PEO chains 
experience in ethanol. 

Fig. shows the composition differences Sci — c^a — 
c^b between chambers A and B divided by the tempera- 
ture difference ST for the three components (PEO, water, 
ethanol) as a function of water content of the solution. 
For PEO in a single solvent, ethanol or water, the com- 
position difference 8c p of PEO is balanced by the compo- 
sition difference of the solvent. In general, the three com- 
position differences add up to zero, Sc p + Sc w + Sc s = 0. 
The composition of the solvent is generally different in 
the hot and cold chamber. For a solution with mass 
fraction of water c w ~ 0.49, for example, there is no dif- 
ference in ethanol content of the solution of the warm 
and cold chamber, while water is enriched in the cold 
chamber. Considering the ethanol content of the sol- 
vent, c s = c s /(c s + c w ), in each of the chambers, we find 
in our calculations that the solvent in the warm cham- 
ber is always slightly richer in ethanol than the solvent 
in the cold chamber. Furthermore, there is a range of 
water contents, approximately 10% to 35% by weight, 
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FIG. 8: Soret coefficient of PEO in mixtures of ethanol and 
water at atmospheric pressure as a function of temperature 
for several solvent compositions. The PEO concentration is 
5g/L at 293 K. 



where the difference in solvent composition c Si A — c Sj b 
between the two chambers is approximately constant at 
(c s ,a-c s , b )/^T = 2.5x10- 4 K- 1 . 

Our results for solvent composition reflect the ap- 
proximations in our calculations, in particular the ran- 
dom mixing approximation for solvent-solvent interac- 
tions. When we perform calculations on ethanol-water 
mixtures without polymer, we find a Soret coefficient 
for water that is positive for all compositions and de- 
creases monotonously from 0.0069 K _1 for ethanol-rich 
mixtures to 0.0054 K _1 for water rich mixtures. Exper- 
imental datai^i^*^, on the other hand, show the Soret 
coefficient of water to change sign with increasing water 
content from around 0.006 KT 1 at 50 wt.% to -0.008 K" 1 
at 95 wt.% 

In Fig. [S] we present results for the temperature de- 
pendence of the Soret coefficient of PEO for several sol- 
vent compositions. For high water concentration, the 
Soret coefficient is positive at room temperature, de- 
creases with increasing temperature and changes sign at 
a temperature that depends on the water content of the 
solution. For ethanol-rich solutions, the Soret coefficient 
is negative at room temperature and increases with in- 
creasing temperature. At the highest ethanol concen- 
trations, an upturn of the St versus temperature curves 
is visible. This signals the onset of evaporation, which 
occurs at increasingly higher temperatures for solutions 
with larger water content. A comparison of Figs. El and 
Fig. [8] with the chain-dimension graphs Fig. 0] (a) and 
(b) shows a correlation between solvent quality and ther- 
modiffusion. In general, as the solvent quality increases, 
indicated by an increase in chain dimensions, the Soret 
coefficient becomes more positive. 
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V. DISCUSSION 

In this work, we have presented a two-chamber lattice- 
model approach to determine Soret coefficients of poly- 
mer solutions. A dilute polymer solution is represented 
by a simple cubic lattice occupied by a single polymer 
chain, solvent particles and, in the case of a compressible 
solvent, voids. Exact enumeration results for an isolated 
chain allow us to construct partition functions for the 
polymer-solvent system without invoking a random mix- 
ing approximation for contacts with polymer sites. Inter- 
actions between solvent particles, on the other hand, are 
evaluated in a random-mixing approximation. Enumera- 
tion results are also used to calculate the radius of gyra- 
tion of the chain in solution, which allows us to monitor 
the solvent quality of the solution. In order to investigate 
thcrmodiffusion, we assume that the lattice is divided 
into two non-interacting sublattices of equal size that are 
maintained at slighty different temperatures. For a given 
occupation of the sublattices, the partition function of 
the combined system is a product of the canonical par- 
tition functions of the sublattices. We consider all pos- 
sible distributions of polymer chain and solvent particles 
among the sublattices. The sum of states of the system is 
calculated by adding up the total partition functions for 
all distributions, while average quantities are calculated 
by performing the appropriate weighted sums. The Soret 
coefficient of the polymer is determined from the differ- 
ence between the average polymer concentration in the 
warm and cold chambers. As in the earlier work on heat 
of transport (cf. Denbigh^) kinetic energy contributions 
are neglected in our calculations. However, we do not 
approximate the heat of transport by a difference in po- 
tential energy. Instead, the probability to find the poly- 
mer in the warmer of the two chambers can be related to 
a difference in average internal energy that reflects both 
enthalpic and entropic contributions. 

We considered first a lattice model for a dilute solution 
of a polymer in an incompressible, single-component sol- 
vent. In this case, the canonical partition function of a 
short chain (N c — 17) in solution and the sum of states 
for the two-chamber system are calculated exactly. We 
derived a simple relationship between the probability to 
find the polymer in the warmer chamber and the internal- 
energy difference U nop ~ U po i, where U nop represents the 
internal energy of a solvent filled lattice and U po \ is the 
internal energy of a lattice containing both solvent and 
polymer. The resulting expression for the Soret coeffi- 
cient shows how our approach relates to earlier work us- 
ing the heat of transport concept. It also illustrates that 
it is important to go beyond a random mixing approxima- 
tion when thermal diffusion in polymer solutions is dis- 
cussed. In agreement with experimental results for dilute 
solutions^ the calculated Soret coefficients are indepen- 
dent of concentration. Unfortunately, the short chains 
considered in this work do not allow us to investigate the 
experimentally observed^ scaling behavior of St with 
molecular mass. We are currently working on a simula- 



tion approach that allows us to treat longer chainaM and 
will investigate the scaling laws in the near future. Our 
results presented in Fig. show that both positive and 
negative Soret coefficients are realized in the incompress- 
ible case and that the sign of St need not be the same for 
all temperatures. A comparison with Fig.^reveals a cor- 
relation between Soret coefficients and chain dimensions 
for two of the three cases presented; the probability to 
find the polymer in the cooler chamber increases as the 
solvent quality increases. This trend appears to be typi- 
cal and, as discussed below, may provide an explanation 
why negative Soret coefficients for polymer solutions are 
so rarely observed. 

In order to describe PEO chains in a mixed solvent of 
ethanol and water, we extended the simple lattice model 
by introducing two types of solvent particles and adding 
voids to account for compressibility. Hydrogen bonding 
between PEO and water molecules is modeled through 
an orientational degree of freedom of the water particles. 
One of the six faces of an elementary cube represent- 
ing a water particle has strongly attractive interactions 
with PEO sites, while the PEO interactions of the re- 
maining five faces are the same as those between ethanol 
and PEO. The lattice model for dilute solutions of PEO 
in ethanol/ water mixtures has eleven system-dependent 
parameters that are determined from static properties 
as described in Appendix 1X1 and summarized in Table [I] 
With these parameters, the model reproduces some of the 
important thermodynamic properties of the system. In 
particular, for PEO in mixed solvents, the solvent quality 
as monitored by the radius of gyration increases as the 
water content of the solution increases^Siii. Similarly, in- 
creasing the temperature increases the solvent quality for 
mixtures with low water content. For mixtures with high 
water content, on the other hand, increasing the tem- 
perature reduces the solvent quality in agreement with 
observations on PEO in water (cf. Ref. l3Cf> . 

Our two-chamber approach allows us to calculate Soret 
coefficients of PEO for given temperature, pressure and 
composition of the solvent. As expected for dilute 
solutions^ the results are independent of the polymer 
concentration. In qualitative agreement with experimen- 
tal data of De Gans et aZii2*ii and Kita et alml, the calcu- 
lated Soret coefficients are negative for solutions with low 
water content and positive for solutions with high water 
content. While we do not expect quantitative agreement 
between our simple lattice model calculations and exper- 
imental data, we would like to discuss some future work 
that may help us gain a better understanding of ther- 
modiffusion in this system. The evaluation of Soret coef- 
ficients with the aid of Eq. is valid only when pressure 
differences between the chambers can be neglected. We 
observe very small differences between the pressures in 
the two chambers and do not expect them to contribute 
significantly, but we are investigating means to correct 
for this effect. Interactions between different polymer 
chains are not considered in our calculations, which ap- 
ply to dilute solutions. However, the polymer solutions 
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used in the experimentsiSiii^ were probably in the semi- 
dilute regime. In typical semi-dilute solutions, the mag- 
nitude of the Soret coefficient decreases with increasing 
concentration^ It will be interesting to investigate the 
effect of chain-chain interactions on negative thermodiffu- 
sion coefficients. Significant thermal diffusion of the sol- 
vent particles is an issue that complicates the comparison 
between theory and experiment for intermediate water 
concentrations^ Eq. J7J for the Soret effect in ternary 
mixtures suggests that several diffusion coefficients need 
to be considered carefully for an experimental determi- 
nation of the Soret coefficient in this regime. Hydrogen 
bonding between solvent molecules may also affect the 
thermal diffusion of the polymer. Our current calcula- 
tions account for specific interactions only between PEO 
and water molecules. We will extend our model to ac- 
count for hydrogen bonds between solvent molecules and 
investigate how it affects the Soret coefficient of the poly- 
mer. 

Finally, we have also investigated the temperature de- 
pendence of the Soret coefficient of PEO and found that 
a change in temperature may induce a change in sign of 
St. A comparison of the results for the radius of gyration 
of the chain and the values of the Soret coefficients re- 
veals a close relationship between the solvent quality and 
the partitioning of the polymer between the chambers. 
In general, the Soret coefficient of the polymer becomes 
more positive as the solvent quality increases. A typical 
experiment on polymers in good solvents is thus expected 
to yield positive Soret coefficients. We expect negative 
Soret coefficients to be observed for polymers that would 
be insoluble were it not for specific interactions between 
solvent molecules and sites on the polymer. It appears 
that both polymer systems for which negative Soret co- 
efficients have been observed, the solutions of PEO in a 
mixed ethanol/ water solventifliiiiS and the solution of 
poly(vinyl alcohol) in water belong to this category. 
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APPENDIX A: DETERMINATION OF 
SYSTEM-DEPENDENT PARAMETERS FOR 
PEO IN ETHANOL WATER MIXTURES 

Our lattice model for PEO in ethanol/water mixtures 
has eleven system-dependent parameters. Seven of these 
parameters describe the thermodynamics of the pure 



components while the remaining four are needed for the 
description of the mixtures. 

For each of the pure components (ethanol, water, 
PEO) we estimate interaction energies en and the vol- 
ume Do per lattice site from a comparison of tabulated 
values^Si for the density and the thermal expansion co- 
efficient at standard temperature and pressure with val- 
ues calculated from lattice-fluid equations of state^^SiS, 
which treat nearest neighbor contacts in random mix- 
ing approximation. For water and ethanol we employ 
a simple lattice fluid model, the Sanchez and Lacombe 
equation-of-state applied to monomers 



P 



1 



-zee 



RT 



ln(l 



(Al) 



where R is the ideal gas constant and vq is the volume 
per mole of lattice sites. 4> is the filling fraction of the 
lattice and related in this work to the mass density p by 



svqP 
M ' 



(A2) 



where M is the molecular mass of the component and s is 
a scale factor that is discussed below. Since the tabulated 
values for the properties of PEO^I correspond to large 
molecular mass, we use an equation of state by Costas 
and Sanctuary^ in the limit of infinite chain length 



P 



2v \ z - 24> J v 
| zRT to / _ 2 
2v \ z 



(A3) 



Comparing calculated and tabulated values of thermo- 
dynamic properties for each substance individually, we 
arrive at three different values for the volume vq. For 
a description on a common lattice, we chose a value for 
v and scale factors Si that allow us to reproduce the 
tabulated densities at standard conditions and obtain 
reasonable estimates for the thermal expansion coeffi- 
cients. The resulting values are presented in Table [I] 
The scale factor for ethanol is chosen as s s = 1 and im- 
plies that each ethanol molecule occupies one lattice site. 
The scale factors for water and PEO are smaller than 
unity. For water, there are l/s w w 3 molecules associ- 
ated with each site. For PEO, our short chains (N c = 17) 
correspond to physical chains with N c /s p ~ 27 repeat 
units. Since each repeat unit has a molecular mass of 
M mono = 44.1 g/mol, the molecular mass of our chains 
is M p = 7V c M mono /s p sa 1187 g/mol. 

The only additional parameters required for the de- 
scription of mixtures of the three substances are the 
mixed interaction energies. Since we are using a very 
simple model to describe the solvent, we employ Bcrth- 
elot's geometric mean combining rule e ws = — -y/e W wess 
for interactions between ethanol and water. This ap- 
proximation leads to deviations between calculated and 
tabulated^ values for the density of less than 2% for all 
compositions at a temperature of 293 K and a pressure 
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TABLE I: System-dependent parameters for the 
PEO/ethanol/ water system 

lattice site volume vo = 5.255 x 10 _5 m 3 /mol 





ey in J/mol 


scale factors 


ethanol 


e ss = -2306 


s s = 1 


water 


e ww = —3306 


s w = 0.3362 


PEO 


e PP = —1153 


s p = 0.6318 




mixed interactions 




ethanol/water 






PEO/water 


fpw.n = 2660 






e P w,s = -8020 




PEO/ethanol 


e ps = 2660 





TABLE II: Comparison of tabulated and calculated val- 
ues for the density p and thermal expansion coefficient a at 
a temperature of 293 K and a pressure of 0.1 MPa 





P (kg/m 3 ) 


a (K- 1 ) 




tabulated 


calculated 


tabulated 


calculated 


ethanol 


789 


788.8 


1.43 x 10~ 3 


1.79 x 10~ 3 


water 


998.2 


998.5 


2.10 x 10~ 4 


3.24 x 10~ 4 


PEO 


1125 


1125.0 


7.2 x 10" 4 


3.83 x 10" 4 



that the chains are not collapsed at room temperature, 
indicating that the short chains are soluble in ethanol. 
However, as in the case of the long chains, solvent qual- 
ity improves with increasing temperature over a large 
temperature range. 

The interactions between PEO and water are described 
by two interaction parameters. Following Jeppesen and 
Kremer— the non-specific interactions e pw;n are assumed 
to correspond to bad-solvent conditions and are approx- 
imated here by the value for the ethanol-PEO interac- 
tions. The specific (hydrogen-bonding) interactions e pw;s 
are strongly attractive. As the remaining free param- 
eter, they determine the solubility of PEO in mixtures 
of ethanol and water. Experiments suggest that the 
transition between poor and good solvent conditions oc- 
curs for mass fractions of water between 5% and 10%. 
We estimate a value for e pw;s from the condition that the 
calculated value of in a mixture with c w = 0.1 is equal 
to R 2 J6*) = 4.13. 



APPENDIX B: DETERMINATION OF LATTICE 
SIZE AND OCCUPATION NUMBERS 



of 0.1 MPa. A better description of ethanol and wa- 
ter can be achieved when specific interactions between 
the molecules are included. We are currently working on 
such an extension of our model. 

With e ss and e pp fixed, the mixed interaction energy e ps 
determines the solubility of the polymer in ethanol. The 
miscibility of polymers is strongly dependent on chain 
length. Unfortunately, solubility data for the short chains 
of our calculations (M p sa 1187 g/mol ) are not available. 
The long PEO chains (M w 236, 000 g/mol) employed in 
the experiments^ are insoluble in ethanol at room tem- 
perature and become soluble only for temperatures above 
323 K. We obtained an estimate for e ps for long PEO 
chains with the aid of the Born-Green- Yvon lattice model 
of Lipson and coworkers^2iiii22. Employing the pure 
component parameters determined above, we adjusted 
e ps until the calculated phase transition temperature was 
near 323 K. The resulting value is e ps = —1141 J/mol. 
This value makes ethanol an excellent solvent for the 
short chains, in agreement with experimental observa- 
tions, but, unfortunately, it does not allow us to use our 
short-chain calculation to explore the effect of changing 
solvent quality with temperature and solvent composi- 
tion. In order to find a more suitable interaction param- 
eter, we use the radius of gyration squared, i?g, as an 
indicator for solvent quality. We estimate a value of e ps 
from the condition that the calculated value of R 2 at a 
temperature of 323 K is equal to R^(6*) — 4.13, the chain 
dimensions of an isolated chain at the 9 temperature of 
the infinite chain. This value for e ps is included in Tabled 
and used in our calculations. Fig. 0] (b) includes a graph 
of calculated chain dimensions as a function of tempera- 
ture for PEO in ethanol. A comparison with Fig.^shows 



In order to calculate properties of PEO solutions at 
given temperature, pressure, and composition, the cor- 
responding lattice size N and occupation numbers JVj, 
i € {p, s,w} have to be determined. In the results pre- 
sented here, we start at room temperature, T = 293 K, 
and atmospheric pressure, P ps 0.1 MPa, and specify the 
PEO concentration as mass per unit volume. This sets 
the lattice size, since the lattice contains a single polymer 
chain occupying 7V C = 17 lattice sites 

N=^/ Pp , (Bl) 

where M p = N c M mono / s p sa 1187 g/mol is the molecu- 
lar mass of the chains and p p is the PEO concentration 
in g/L. If two solvents are present, we specify the oc- 
cupation number for one them, say water (N w ). The 
remaining N — N c — N w lattice sites are then divided be- 
tween ethanol (N B ) and voids (N v ) so that the pressure 
calculated according to Eq. (|2(j[l is approximately equal 
to atmospheric pressure. Since the lattice occupation 
numbers are integers, the targeted pressure value cannot 
be reached exactly. We use an iterative algorithm that 
identifies occupation numbers N s and N v corresponding 
to the closest calculated pressure above the target pres- 
sure. The uncertainty introduced into our results by the 
integer nature of the occupation numbers decreases as 
the size of the lattice increases. Since the number of lat- 
tice sites is set by the PEO concentration, and since our 
results for the Soret coefficient of PEO are independent 
of its concentration, large lattices can be used to reduce 
the uncertainty in the calculated values. Unfortunately, 
the computation time grows rapidly with increasing lat- 
tice size, so that a compromise has to be found. In this 
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work, we use p p = 5 g/L at T = 293 K and target pres- 
sure P = 0.1 MPa. This corresponds to a lattice size of 
N — 4516 sites and yields calculated pressure values from 
0.13 MPa for ethanol rich mixtures to 0.4 MPa for water 
rich mixtures. In order to perform calculations at ele- 
vated temperature or pressure, we retain the occupation 
numbers for the filled sites, 7V C , N w , and N s , but change 
the total number of lattice sites by adding or removing 
voids to reach the desired conditions. 

The concentration of each of the components in the 
system is related to the lattice occupation through 



Nv ' 



Pi=P, ie{s,w,p} 



(B2) 



where 



N s Ma 



N r M,r 



(B3) 



and M s = 46.07 g/mol and M w = 18.0 g/mol are the 
molecular masses of ethanol and water, respectively. The 
mass fractions a of the components are given by 



Pi/P- 



(B4) 
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